Perturbative generation of ^13 from tribimaximal neutrino mixing 
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Abstract 



Solar and atmospheric neutrino oscillations are consistent with a tribimaximal form of the mixing 
matrix U of the lepton sector. Exact tribimaximal mixing leads to ^13 = 0. Recent results from the 
J> . Daya Bay and RENO experiments have established a non-zero value of ^13. Keeping the leading be- 
I haviour of U as tribimaximal we perform a model-independent perturbative calculation to incorporate 
\^ ' a non- vanishing ^13. We identify the nature of the perturbation matrix and consider the possibility of 
■ the solar neutrino splitting also resulting from it. We calculate up to first order in perturbation theory 
Tij- I and evaluate the deviations proportional to sin0i3 while including CP-nonconservation. Finally, we 
' briefly discuss a gauge model where such an addition to the neutrino mass matrix arises through 
one-loop effects. 
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Experimental data of solar, atmospheric, accelerator, and reactor neutrinos [T] translate to information 
about neutrino masses and mixing which can be summarised as [21 [3]: 
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(7.59 ± 0.20) X 10"^ eV^ 6*12 = (34.4 ± 1.0)°, 
(2.46 ± 0.12) X 10"^ eV^ 623 = (42.8l|:^)° 
(5.6t|?)°, (^unknown. 



(1) 



These values of the mixing angles are consistent with a mixing matrix of tribimaximal form [3] , 
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which predicts the third mixing angle ^13 to be exactly vanishing. 

Of late, the situation has taken a different turn. Results from the Double Chooz collaboration and 
more recently the Daya Bay [6j experiment indicate that ^13 is, in fact, inconsistent with zerc0 by 
more than 5a - sin^ 2^13 = 0.092 it 0.016 (stat)± 0.005 (syst) [6]. Therefore, in view of these significant 
findings, it has to be concluded that the simple-minded tribimaximal picture fails to adequately 
capture the observed neutrino mixing. The smallness of ^13 compared to the other two mixing angles 
encourages us to examine here whether the former could arise from a small perturbation on the basic 
tribimaximal structure and could lead to a realistic neutrino mixing matrix. 

We work in a flavor basis in which the charged lepton mass matrix is diagonal!. If the left-handed neu- 
trino Majorana masses be mi, 1712, then from eq. ([2]) the mass matrix M^, satisfying tribimaximal 
mixing, when expressed in the flavor basis has the general form. 
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(3) 



where we have set 



mo = {mi +m2 + m3)/3, A32 = {m^ - 7712) and A31 = (7713 - mi). (4) 

Ah initio, the mass eigenvalues, mi, 7712,7713, can be complex in which case they can be rendered real 
and positive by a diagonal phase transformation, D = diag{e^^^ ,e^^'^ ,1), where the Aj are Majorana 
phases, which do not affect neutrino oscillations. 

We approximate A32 ~ A31 = A, which is not unreasonable since IA32I 3> A21 = (m2 — mi). A sets 
the scale for atmospheric neutrino oscillation^. We start with this limit and write the unperturbed 

^Very recently the RENO collaboration has measured sin^ 26*13 = 0.113 ± 0.013 (stat) ± 0.019 (syst) 7 . 
^This fixes the singlet right-handed charged leptons and the left-handed lepton doublets in flavor space. In this basis, 
mixing in the lepton sector is determined entirely by the neutrino mass matrix. 
•^A is positive (negative) for the normal (inverted) ordering of neutrino masses. 
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mass matrix in the flavor basis as: 






( " mo + ^ # I . (5) 

T mo + f , 

At this level, m^^ = m!^^ = mo — ^ and m'^^ = niQ + ^ and the solar mass splitting is absent. Our 
goal is to also generate this splitting through the same perturbation hamiltonian that is responsible 
for 7^ 0. We take m'^\rn~2 \ ^'^d m!'^^ to be real and positive. 

The purpose of this paper is not to explain how emerges from a fundamental model; even though 
there is no doubt that we consider it as the dominant part of the neutrino mass matrix. There are many 
models from which one can obtain the tribimaximal form of the mixing matrix Our discussion 
below will be independent of the specific mechanism by which arises. 

In terms of the three mixing angles and the complex phase 5 the Pontecorvo, Maki, Nakag awa, Sakata 
(PMNS) mixing matrix is conventionally parametrized as, 

C12C13 S12C13 
U = I -S12C23 - ci2S23Si3e*'^ C12C23 - si2S23Si3e*'^ S23C13 | • (6) 

S12S23 - Cl2C23Sl3e*'^ -C12S23 - Sl2C23'Sl3e*'^ 




As noted, the tribimaximal mixing matrix, in eq. ([2]), fixes the element U^^ = 0. The role of a 
non- vanishing C/e3, or equivalently ^13, is manifold. It is essential for CP-nonconservation in neutrino 
oscillations and may be invoked to explain leptogenesi^. Also, ^13 7^ will be similar to the quark 
sector where mixing between all three generations and CP- violation is a well- verified result, though 
the mixing angles in the two sectors are vastly different. For CP-violation, of course, both ^13 and 
the complex phase 6 should be non-vanishing. Besides, a reasonably large ^13 opens the door for an 
easier measurement of the neutrino mass ordering, i.e., the sign of Amg^. 

A large number of attempts have been made to generate ^13 7^ in diverse ways starting from an 
initial tribimaximal form. Some of these are the following. A perturbative analysis in which one of 
the columns or rows of is left unchanged has been examined in [9]. An alternative which involves a 
sequential 'integrating out' of heavy neutrino states has been proposed in |10] . Another approach has 
been to parametrize the deviation from the tribimaximal form in a particular way Deviations 
from tribimaximal mixing due to charged lepton effects and Renormalization Group running have 
been other directions of study [12J. Alternative explorations have been based on the A{4) symmetry 
in |13t |T3] , and on other discrete symmetries in |15t 116] . 

Our strategy here is to use perturbation theory to identify the structure of the Majorana mass matrix, 
M = Af° + M', where M' < M^, so that 6*13 and the solar mass splitting are obtained. Both and 
M' will be symmetric and could, in general, be complex. However, as obtained in eq. ([5]) from 
the tribimaximal mixing form is real and symmetric, i.e., hermitian. We will consider the cases of real 
and complex M' separately. 

''The Majorana phases alluded to earlier could produce CP-violation in AL — 2 processes. 
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For our later discussions the eigenstates of M^, the unperturbed mass eigenstates, in the mass basis 
are found useful. These are simply: 
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(7) 



of which the first two are degenerate. So, the basis vectors tp^^ and ip^^ are not unique and are chosen 
with the knowledge that they reproduce the correct solar mixing. The physical basis is fixed by the 
perturbation. When we discuss lifting of the degeneracy, we consider M' to be such that ip^^ and tp^^ 



are its non-degenerate eigenstates: < ipf'^lM'l'tpj 
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m- 0; 



= 1,2), with m^^ ^ rri^' ■ We also 
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take (M')33 = in this mass basis, so what remain are (M')i3 and (M')23 to which we will first turn. 

It is helpful to bear in mind that eigenstates in eq. ([7]) when expressed in the flavor basis are simply 
the columns of , eq. ([2]), namely. 
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The goal we have set ourselves is to obtain as the perturbed mass eigenstates, when written in the 
flavour basis, the columns of the matrix in eq. ([6]) with ^13 7^ 0. To this end, initially, let us take M', 
which is symmetric, to be real and therefore hermitian. Needless to say, this may generate a non-zero 
^13 but will have no CP-violation and hence yielcH 5 = 0. For the perturbation expansion we retain 
terms up to linear in S13. To flrst order we have. 
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The coefficients O^j are real in this case. In the mass basis O 



(i / 3). 

is proportional to Mij. 
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The eigenstate, '03 j should correspond to the third column of the mixing matrix U in eq. ([6|) with 
6 = 0. O31 and O32 are readily determined using eq. ([9]) in the flavor basis. Written explicitly we get 
the matrix equation, 
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3S13, 



By inverting the above equation one obtains, to order linear in S13, O31 = y |si3 and O32 = 
where maximality of the atmospheric mixing angle (§23 = C23 = l/\/2) has been used. This translates 
to M^g = w IsiaA and M23 = w ^si3A in the mass basis. 



Note that a negative S13 with 5 = is equivalent to a positive S13 and 5 — n. 
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To extend this discussion to the case of 5 7^ 0, we have to bear in mind that now M' is complex 
symmetric and not hermitian and the same holds for the total Majorana mass matrix M = + M' . 
The columns of the mixing matrix U (eq. ^) are eigenvectors of M+M = M^^ + M^^ M' + M'^ , 
where we have dropped a term which is 0(M')^. To proceed, we recall that is hermitian and 
therefore the eigenstates of the unperturbed M^"^ are the same il^f'^ considered earlier - eq. ([8]) - 
but now corresponding to eigenvalues (m^'^^)^, and {m'^^)'^ . In place of eq. (jlOp we have 

<V.f |(MOtM- + M-tM°)|^f > 
= r^^^^^JJ^^ = / 3) , (12) 



mg j -ym. 



which is to be used in eq. ([9]) now. Requiring that be reproduced to first order and using the 
appropriate variant of eq. (fTTI) we get in this case O31 = ^J^si^e~'^^ and O32 = y'^sise"*''. 

To relate the above to the elements of the perturbation M' one notes: 

< KMOtM' + M'tMO)|V'f ^ >= mf < ^j°)|M'|V'f ^ > + mf^ < i;f\M'^\i;f^ > , (13) 

and thus in the mass basis 

2 / inWi' 



O 



3j 



(mf) - (mf) =mf\M'),;+mf\M%, {j ^ 3) , (14) 



where the symmetric nature of M' has been used. Writing = y ^s'^^A e^"^ and Af23 = y I'^'is^ c*"'^ 
using eq. (114^ one finds to leading order in A/mo 

(5 = tan^^ f tan(/) ) , and sis = /(</))s'io, (15) 

V2mo y 

where 

'(mf )2 + (mf )2 + 2mf mf cos 2 J 
/('^) = ^ -W. u^. — ■ (16) 



/ (0) , (0)^ 

[rrii + mg ) 

The approximate formulae in (llSp indicate that S13 < with the equality holding only when </> = 0, 
and though the range of - which is {0, 27r} - is the same as that of 6 the latter is suppressed 
compared to the corresponding (j). The suppression is higher as the neutrino masses approach the 
quasi-degenerate regime (A ^ mo). 

So far we have concentrated on obtaining ^13 7^ through a perturbation starting from the tribi- 
maximal form. Now we consider the solar mass splitting. We choose the perturbation such 
that {M')i2 = {M')2i = 0. The first order corrections to the neutrino mass are obtained from 
mf = < ipf^\M'\ipf''^ >. We demand that the following mass corrections arise at this order 

mf = mf = and mf 7^ 0. (17) 

In the mass basis this implies that out of the diagonal elements only {M')22 7^ 0. Such a correction 
ensures that a nonzero solar mass splitting 1112 — mi = mf is induced. Solar neutrino observations 
establish Am^^ = (^-2)^ — (mi)^ is positive. 
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Putting all this together we have for the full perturbation matrix in the mass basis as 



M' = s'lgA 



2 i4 

3 ^ 
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The dimensionless parameter x is fixed by the solar splitting. In general it can be complex implying 
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(20) 



A is a Majorana phase of V2 which arises from the perturbation. 

There are thus two real parameters introduced here: |?7i2^^| and x- For any phase angle x demanding 
that the solar splitting is correctly obtained determines l^n-g^^l provided rni^ is known, i.e., the mass 
ordering is specified and the mass of the lightest neutrino, m, is given. Thus x-, ■s'13 and suffice to 
fix the full perturbation matrix M' . 

Using M' in eq. (|18p and degenerate perturbation theory we get for the mixing matrix with 5 7^ 0: 
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I Size'' 
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(21) 



Us^o is consistent with the observed mixing angles and is unitary up to order S13. The non-zero 
CP-phase 5 brings the lepton sector in line with the quarks, where CP-violation has been established 
for long. 5 is usually invoked for processes such as leptogenesis. A matrix of exactly the form of Us^o 
has been discussed in [18] from a different motivation and its consistency with the experimentally 
required mixing angles noted. 

The basis independent measure of CP-violation, the leptonic Jarlskog [19] invariant, arising from Us-^q 
(eq. (EU) is 



J = lm[UeiU^2U:2U. 



fill 



3a/2 



si3 sin 5 



3V2 



s'l3f{<P)- 



SUK 



COS^ 
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V 2mo ) 



(22) 



sm 



signifying that both and (j) have to be non-vanishing in order for CP-violation to be present in the 
lepton sector. Moreover, in the quasi-degenerate regime the observation of CP-violation is less likely. 

The above discussion is valid when the solar mass splitting and the mixing angle ^13 are unrelated. In 
the following we do not examine mass matrices of the associated general form - eq. (jl8p . Nonetheless, 
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we make one passing remark. It would not be unreasonable to expect that the different non-zero terms 
of the perturbation matrix (jlSp are rougly of similar order. We may then expect x ~ 0(1). Recalling 



eq. p9|) one has the order of magnitude estimat^ S13 ~ 0[{Am2i/ A.m'^j^){m 



(0) 



, (0)x // (0) , (0) 
+ m\ )/(m2 +m\ 



The measured values of Am'^i and |Am3^| are known. We illustrate two extreme limits 
C^) ^^^\mi!^^ implies S13 



)]■ 
normal 



ordering with m 



o\w- 



(0) 



/2mf'^)] while for the inverted ordering 



3 » m\ ' ,1712 ' implies S13 ~ ^[iU "(m-^ 
with m^^^ <C 771^^^777-2^^ ^13 ~ ©[lO^^]. This is the general expectation if both ^13 and the 

solar mass splitting arise from the same perturbation of the tribimaximal mass matrix. 



We now identify a special limit when the perturbation mass matrix is of a texture which can be realised 
from a simple model and where S13 gets related to Am^i resulting in restrictive predictions. To relate 
to mass models it is more convenient to first rewrite M' in the flavor basis. We find from eq. (llSp 
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1 i4 
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in flavor basis. 



(23) 



Here the first matrix on the right-hand-side is responsible for ^13 and the second for ATn^j^. 

We see from eq. (|23p that, aside from the diagonal part which is proportional to the identity matri^ll 
and can be subsumed in M", the perturbation is of the form: 



M' 




in flavor basis, 



(24) 



where A, B, and C are complex in general. Such a texture of M' can follow from a Zee-type model 
pO] as we discuss later. In such models (M')^^ is proportional to (777^ — 777 



of the charged lepton a. As rrij- ^ m 



where m„ is the mass 



Hi 



777p. 



unless other couplings are of vastly different order from 
each other, one must have ~ C S> A. Such a form of the mass matrix can be reproduced by the 
choice -^^^^ + X = e, where e is small, when the perturbation matrix eq. ()23p reduces to: 



M' 





(25) 



The special case e = is quite predictive. From eq. (|T9|) this requires 



V2 



'13 



A e" 



m. 



(1) 



(26) 



Thus 5^3 = V2 V3|A| and 

= TT + X {(p = x) for the normal (inverted) ordering. 

Due to these relationships, x ^-nd besides determining |7n2^^| now also fix S13 and 5 through eq. 
()15p . As noted, mf'^ is known when the mass ordering and the lightest neutrino mass, 777, are fixed. 
So, for any mass ordering the two remaining parameters are x ^-nd the lightest neutrino mass. 



This result is only indicative. The full flexibility of variation of <^ and x is not taken into account. 
'^Such a piece proportional to the identity does not affect the mixing and makes a constant contribution to all three 
neutrino mass eigenvalues. 
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' r 1+ > 

— ► ^ < X — ► ^ — 



f 

Figure 1: One loop contributions to the neutrino Majorana mass matrix M' in tlie Zee model. For M'^^ the 
dominant contribution is proportional to while for M'^^ and M'^^ it is proportional to m,-- 



We show now that x ^^id rh can be chosen such that one has consistency with both the solar mass 
splitting and the measured ^13. In our discussion below we use the central values of the atmospheric 
and solar mass splittings from eq. (JI]) and seek an acceptable ^13. We do not attempt an exhaustive 
listing of the entire consistent ranges of the parameters in this work but rather present some typical 
solutions for both mass orderings. 

We find that with e = 0, in the normal mass ordering case {m'^^ = m) taking mf'^ = 10~^ eV and 
115° < X < 137°, sin2 26li3 varies from 0.057 to 0.130, which includes the la experimentally allowed 
range, while the Jarlskog CP-violation parameter J remains more or less constant around -0.026. The 
replacement x ^ (^tt — x) with m^i^ fixed keeps sin^ 2^13 unchanged and replaces J by — J. 

For the inverted mass ordering (m^'^^ = m — A), on the other hand, taking < m.g'^^ < 10~^ eV and 
X ~ 94° one obtains sin^ 26I13 ~ 0.054, which is allowed at 2a, with J ~ -0.028. For mP = 3 x 10"^ 
eV and x ~ 101° o^is has sin^ 2^13 ~ 0.09 within la and J ~ —0.029. For any , replacing x by 
(27r — x) results in the same sin^ 2^13 but J changes sign. 

We now briefly note how M' of the texture in eq. (j24p can follow from a Zee-type models]. It bears 
repetition that here the Zee model provides a subleading contribution, M' , to a leading tribimaximal 
mass matrix of a different origir^j. The Zee model has a simple SU{2)l x U{1)y invariant structure. 
For this, a second scalar SU{2)l doublet and a charged singlet scalar x'^ are introduced. The latter 
couples to a pair of lepton doublets, where the coupling fap is antisymmetric in the generation index. 
Likewise due to SU (2) antisymmetry the charged scalar also couples to a pair of Higgs doublets hu 
and hd antisymmetrically. In this model a contribution to the neutrino mass ~ M' ~ arises radiatively 
from one loop diagrams such as Fig. [T]and can be expressed as: 

</3 = liTl/^K - • (27) 



'In this case, rh = nn^ and A < 0. 

'An alternative way to generate a non-zero 6^13 using the Zee model has been examined in [21], 

'Models can be constructed which accommodate both M° and M' . An A(4) based example can be found in [14]. 
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Here fajS is the antisymmetric coupling in fapLaLf^x^ where L is the left handed lepton doublet. Also, 
^ is the trilinear scalar coupling in jihuhdXi a typical scalar mass, and / a dimensionless factor 
arising from the loop integral. VuiVd are the vacuum expectation values of the two Higgs doublets hu 
and hd- The vertex / violates lepton number by two units. This diagram gives rise to a mass matrix 
M' which is off-diagonal and symmetric, as required. {M')i2 can be neglected compared to (M')i3 
and (M')23 because the latter will receive contribution from diagrams with a r- leptor0. Thus, the 
correction obtained in this fashion is naturally of the desired form with e ~ 0. Further, the coupling 
fa^ can be complex which can lead to an M' of the form in eq. (j24p - which is complex symmetric. 
The interference of M' with the matrix - eq. (jl2p - leads to CP-violation in the neutrino sector. 
It is worth bearing in mind that M' is suppressed compared to the leading term in by ©(sisA/m). 
Taking S13 ~ 0.1, A ~ 0.1 eV and ^ ~ 100 GeV, unless other factors in eq. ([27|) are tuned to suppress 
the contribution, one requires Mg ~ 0(10^ GeV), which puts the additional scalars of the model 
beyond the reach of the current experiments. 

In conclusion, we have shown that ^13 consistent with experiments, a CP-phase 5, and the solar mass 
splitting can all be the outcome of a specific perturbation to a basic neutrino mass matrix, the latter 
associated with tribimaximal mixing. This leads to a non-zero Jarlskog invariant and opens the door 
for CP-violation in the lepton sector. In particular, a constrained version of this perturbation relates 
the neutrino Majorana phase to the solar mass splitting as well as ^13 and 5. Some sample solutions 
which meet all requirements have been presented. We have provided an example where the requisite 
perturbation contributions to the neutrino Majorana mass matrix can arise from a Zee-type model 
through radiative corrections. 
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